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Abstract Thermal biology traits reflect thermal 
adaptations to an environment and can be used to 
infer responses to climate warming in animal species. 
Within a widespread genus or species, assessing 
the latitudinal or altitudinal gradient of thermal 
physiological traits is essential to reveal thermal 
adaptations and determine future vulnerability to 
climate warming geographically. We determined the 
thermal biology traits of a cold-climate distributed 
lizard, Takydromus amurensis, and integrated 
published thermal biology traits within the genus 
Takydromus to reveal a preliminary geographical 
pattern in thermal adaptation. The mean selected 
body temperature (cloaca temperature; T,,,,), critical 
thermal maximum (CT,,,,), critical thermal minimum 
(CT 

sprint speed; Taw of T. amurensis were 32.6, 45.1, 3.1, 


and optimal temperature for locomotion (i.e., 


an) 
and 33.4 °C, respectively. The resting metabolic rates 
of T. amurensis were positively related to temperature 
from 18 °C to 38 °C. We compared the traits of 
tropical T. sexlineatus, subtropical T. septentrionalis, 
and T. wolteri with T. amurensis and found that the 
CT jax and thermal tolerance range (the difference 
between CT a and CT,,;,; TTR) increased toward 


min? 
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high latitudes, whereas CT „in increased toward low 


latitudes in these four Takydromus lizards. According 


min 


to this preliminary pattern, we speculate the species at 
medium and low latitudes would be more vulnerable 
to extreme heat events caused by ongoing climate 
warming. We highlight the importance of integrating 
thermal biology traits along geographical clues, and its 
potential contribution to evaluate the vulnerabilities of 
species in the context of climate warming. 


Keywords counter gradient, CT mav CT miw thermal biological 


muy 


trait, thermal tolerance range, T „p Takydromus 


sel 


1. Introduction 


Ongoing climate warming has imposed novel stressors 
on animals (e.g., Thomas et al., 2004), which has affected 
their populations in a number of ways. For example, 
many species have altered their distributions, and some 
populations face risk of collapse (e.g., Pounds et al., 1999; 
Wilson et al., 2005). Recent studies have shown that the 
thermal safety margin (difference between environmental 
temperature and critical maximum temperature, 
CT max) Of low-latitude animals is narrower than that of 
animals living in mid and high latitudes, and the thermal 
biological characteristics of low-latitude animals are 
more susceptible to temperature variation, indicating 
that low-latitude animals are more vulnerable to climate 
warming (e.g., Dillon et al., 2010; Giisewell et al., 2017; 
Sunday et al., 2014). Thus, a detailed understanding of 
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the thermal biology traits and the thermal responses of 
animals is not only useful for studying adaptations to 
thermal environments, but also necessary for determining 
the vulnerability to thermal variation caused by climate 
warming. 

Thermal biology traits have been shown to be 
important proxies of the thermal adaptation of animals 
to their environments, as measured by selected body 
CT nin and the thermal 
dependency of performance such as locomotion (e.g., Huey 
and Hertz, 1984; Huey and Kingsolver, 1989; Angilletta 
et al., 2002a). Unlike endotherms, which can maintain 


temperatures (Tsa), CT navi min? 


a stable body temperature in a wide temperature range 
at great energetic costs (Bennett and Ruben, 1979), the 
body temperatures of reptiles and other ectotherms are 
regulated by thermal interactions with the environment 
and thus are more susceptible to thermal variations and 
their effects upon behavior and physiological processes 
(Huey and Stevenson, 1979; Angilletta et al., 2002a, 
2002b). Accordingly, reptiles have a high risk of extinction 
owing to climate warming (e.g., Sinervo et al., 2010). 
In addition, they are affected by internal and external 
factors (Wang et al., 2017). Through the analysis of the 
thermal biology traits of widespread genera or species, 
inferences can be made regarding the geographical 
patterns and interspecific variation of traits. This can aid 
in determining the vulnerabilities of widespread genera to 
climate change along geographical gradients (Huey et al., 
2009; Sinervo et al., 2010). 

Reptiles can maintain their body temperature within 
a moderate range by thermoregulation, with a range 
of optimal body temperatures for different processes 
(Van Damme et al., 1991; Bauwens et al., 1995). The 
selected body temperatures of a reptile represent the 
body temperatures without biotic or abiotic constraits 
(Huey et al., 1979; Hertz et al., 1993), which can be 
determined by the temperature gradient in the laboratory, 
although extrapolation of laboratory data to field 
populations should be interpreted cautiously. Because 
T,.; represents the “ideal” range of body temperatures 
ectothermic animals aim to maintain through behavioral 
thermoregulation, it is normally different from active 
body temperatures (Licht et al., 1966; Avery, 1978). CT 
and CT 
at which vitality can be sustained, respectively (Hertz et al., 
1983; Gilbert and Miles, 2017). If the body temperature 
is lower than CT,,,, or higher than CT 
lose motor function and suffer cold narcosis (temperature 


max 


are the maximum and minimum temperatures 


min 


reptiles would 


max) 


lower than CT,,;,) or muscle spasms (temperature 
higher than CT,,,,) (Lutterschmidt and Hutchison, 1997; 
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Angilletta, 2009; Camacho and Rusch, 2017). Therefore, 
CT wax and CT 
reptiles in various thermal environments (Hertz et al., 
1983; Deutsch et al., 2008). 

Metabolism is the primary physiological process that 


ee min are used to determine the viability of 


determines an organism’s demand from the environment 
and the allocation of energy among body functions 
(Brown et al., 2004). In ectotherms, the resting metabolic 
rate is often used as a proxy for metabolism and is 
defined as the energy expenditure of animals while in 
the resting, fasting, and nonproductive states (Andrews 
and Pough, 1985; Christian et al., 1999; Berg et al., 2017). 
In reptiles, resting metabolic rates are highly sensitive 
to body temperature, requiring the organism to respond 
quickly to variations in the thermal environment (Ma 
et al., 2018a, 2018b; Sun et al., 2018). As another critical 
function, locomotion is also significantly sensitive to 
body temperature and is related to fitness because of 
its importance for foraging, avoiding predators, and 
choosing mating partners (Bauwens et al., 1995; Chen 
et al., 2003; Sun et al., 2018). In summary, evaluating 
the thermal sensitivity of critical indices (i.e., metabolic 
rates and locomotion) is essential to understand thermal 
adaptation in ectotherms (Angilletta, 2009). 

Takydromus lizards have a large geographical 
distribution, with a significant latitudinal span, ranging 
from cold-temperate to tropical regions. As a species 
that inhabits cold climates, the Heilongjiang grass lizard, 
Takydromus amurensis, is a small lacertid lizard (up to 
70 mm snout-vent length [SVL]) which is distributed 
in northeast China and the boundary of China, Russia, 
and the Korean Peninsula (Xu et al., 2017). Like other 
Takydromus lizards in mainland China, T. amurensis 
inhabits areas at the boundary of forests and grasslands 
or scrub habitats. Females are able to lay one or two 
clutches from May to July, with an average clutch size 
of six eggs (Zhao et al., 1999). T. amurensis is distributed 
in high latitudes, which typically exhibit low average 
temperatures and drastic thermal fluctuations (Figure 
1). Thermal biology traits should be investigated in T. 
amurensis to obtain a comprehensive understanding 
of thermal adaptation in cold-climate ectotherms. The 
thermal biology of several other congener species has been 
previously studied, including T. wolteri, T. septentrionalis, 
and T. sexlineatus, which inhabit subtropical and tropical 
areas, respectively (Ji et al., 1996a; Chen et al., 2003; Zhang 
and Ji, 2004). By assessing the same traits in T. amurensis, 
i.e., those investigated in these previous studies, we can 
improve our understanding of the specific variations 
and preliminary patterns of thermal biology traits in 


Asian 
Herpetological 
Research 


352 


E] Thermal fluctuation 
Average temperature 


Temperature (°C) 


Jan Feb Mar Apr May Jun Jul 


Aug Sep Oct Nov Dec 
Time of year 

Figure 1 Annual temperature environments for Anshan popula- 

tion of T. amurensis, where we collected the lizards. The red line 


indicates the average air temperature along the year, and grey 
area indicates the daily fluctuations of air temperature. 


Takydromus species along latitudinal and altitudinal 
gradients. 
We determined selected body temperatures, thermal 


tolerance (i.e., CT aax and CT nin) and thermally sensitive 


max 
traits including resting metabolic rate and sprint speed 
as focal study traits for the Heilongjiang grass lizard, 
T. amurensis. We aimed to (1) contribute with new data 
on a cold-climate distributed Takydromus species (ie., T. 
amurensis) for the first time; (2) compare thermal biology 
traits among Takydromus lizards with published data 
concerning T. septentrionalis, T. wolteri, and T. sexlineatus; 
and (3) summarize preliminary patterns of variation 
in Tep CT mav CT miy and optimal body temperature for 
performance in Takydromus lizards along a latitudinal or 
altitudinal gradient. 


2. Materials and Methods 


We collected adult T. 
amurensis (n = 13; 4 males and 9 females) in late April 2016, 


2.1. Study species collection 


using either a noose or by hand, from Anshan (see details 
in Figure 1 for annual thermal environments), Liaoning, 
China (41°01 N, 120°7’ E; ~210 m). After collection, we 
transported the lizards back to the laboratory in Beijing. 
The lizards were individually housed in plastic terraria 
(350 mm x 250 mm x 220 mm, length x width x height). 
The terraria were set up with a mixture of moist soil and 
sand, with randomly placed patches of grass. The lizards 
were housed in a temperature-controlled room at 18 
°C under a natural photoperiod, with a supplementary 
heating lamp suspended above one end of the terrarium 
from 06:00 to 20:00. During the heating period, 
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temperatures in the terraria ranged from 18 °C to 40 °C. 
Food (crickets and larval Tenebrio molitor) and water were 
provided ad libitum. After three days of captivity, thermal 
biology traits were started to be determined. 


2.2. Selected body temperatures (T,,,) Selected body 
temperatures (T,,,) were determined according to the 
established protocol, with minor modifications (Shu et 
al., 2010). In brief, T 
terrarium (1000 mm x 500 mm x 300 mm, length x width 


sei Was measured in a custom-made 
x height) placed in a temperature-controlled room at 
18 °C. The terrarium was fitted with grass substrate and 
‘natural’ retreats to mimic the natural environment. One 
heating light bulb (275 W) was suspended above one end 
of the terrarium, creating a thermal gradient from 18 °C 
to 60 °C during the heating period. The heating period 
lasted 14 h daily from 06:00 to 20:00. On the first day, a 
group of lizards (four to five) were introduced into the 
terrarium from the cool end at 16:00 for acclimation 
to the surroundings. Then, on the second day, the body 
temperatures of the lizards were measured twice at 09:00 
and 15:00, with an UNT-325 electronic thermometers 
(UNT T-325, Shanghai, China) by inserting the tip of the 
probe into cloacas. During the measurements, each lizard 
was captured by hand and its body temperature was 
measured immediately (within 30 s), without disturbing 
the other lizards in the terrarium. Then, the average 
of the two body temperature records of each lizard 
was estimated and used as the lizard’s selected body 
temperature. 


2.3. Thermal tolerance The critical thermal minimum 
(CT,,;,) and critical thermal maximum (CT nax) 
temperatures of the lizards were determined in a 
programmed incubator (KB 240, Binder, Germany). We 
cooled (for CT 
determination) lizards from a starting temperature of 


min determination) or heated (for CT,,,,, 


28 °C ata rate of 1 °C per 10 min. During the cooling and 
heating, we continuously monitored the behavior of the 
lizards. When the lizards lost the capacity to respond to 
intense stimulation and could not right themselves after 
being turned over, the body temperatures were recorded 
as CT pin OF CT maw respectively (Zhang and Ji, 2004; Xu 
and Ji, 2006). The lizards were then moved to a 28 °C 
chamber to recover. The protocol stands that if a lizard 


max? 


cannot recover from the cold or heat shock, its record is 
eliminated for further analysis (Li et al., 2017). However, 
in this study all the lizards were recovered from the shock 


2.4. Locomotor performance Locomotor performance 
was estimated using sprint speed. The sprint speed was 
measured at five test temperatures ranging from 18 °C to 
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38 °C (18 °C, 23 °C, 28 °C, 33 °C, and 38 °C), based on the 
reported range of field body temperatures of the species 
(Xu et al., 2017), in a randomized sequence. Locomotor 
performance was measured once per day under one 
temperature treatment (18 °C, 33 °C, 28 °C, 23 °C and 
38 °C), in a sequence order. Before the test, the lizards 
were placed in an incubator at the test temperature for 
approximately 2 h for acclimation. To ensure that body 
temperatures of lizards matched those of the respective 
test treatment, we measured the body temperatures of a 
subset of lizards before starting stimulating them to run 
(Sun et al., 2014). Locomotor performance was tested by 
stimulating the lizard to run through a racetrack, which 
was recorded by an HD video camera (Sony, DCR- 
SR220E, Japan). The racetrack was 1500 mm x 100 mm 
x 150 mm, with intervals marked every 250 mm. Each 
lizard was stimulated by a paintbrush to run twice at each 
temperature with an interval of 1 h for rest. The videos 
were analyzed by Windows Movie Maker. For each lizard, 
the fastest speed through 250 mm for each time was 
recorded, and the average of the two fastest records was 
used as the sprint speed (Sun et al., 2014). 


2.5. Resting metabolic rate (RMR) The resting metabolic 
rates were estimated by respiratory gas exchange rates 
measured at five random test temperatures (18 °C, 
23 °C, 28°C, 33 °C, and 38 °C). Before the test, the 
lizards were fasted for at least 12 h. All lizards were 
acclimated at the correspondent test temperature for 2 
h in an incubator (KB 240, Binder, Germany). Then, the 
lizard was enclosed within the respirometry chamber 
placed in the incubator (KB 240, Binder, Germany). The 
respiratory gas exchange was determined using a closed- 
flow respirometry system with a volume of 281.4 mL 
(Stable System International Inc. Las Vegas, NV, USA), 
and the resting metabolic rate was estimated via the CO, 
production rate using a previously established method 
(Sun et al., 2018). For stabilizing the gas composition of 
the system (i.e., baseline), we opened the system to the 
air for approximately 5 min before changing the system 
to closed-circuit respirometry. Then, the carbon dioxide 
production rates (i.e., VCO,) in the closed circuit were 
continuously recorded for approximately 10 min. To 
minimize the effects of circadian rhythms, measurements 
were conducted from 10:00 to 18:00. The metabolic 
rates were calculated as the CO, production per gram of 
body mass per hour (mL/g/ h), following the equation 
metabolic rate = VCO, x volume/body mass, where VCO, 
is the CO, production rate in percentage (%/h) in the 
closed circuit with a volume of 281.4 mL. 
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2.6. Statistical analysis The normality of distributions 
and homogeneity of variance were tested using the 
Kolmogorov-Smirnov test and Bartlett’s test prior to 
analysis. First, we tested for sex differences in Tep CT maw 
and CT 
as a factor. As there were no significant differences (all 
P > 0.05), we pooled the data of both sexes for all these 
traits. For the analysis of locomotion, we performed the 
Gaussian (R° = 0.473), modified Gaussian (R° = 0.523), 
Lorentzian (R’ = 0.513), Weibull (R’=0.515), and Pseudo- 
Voigt regressions (R° = 0.512) to analyze the sprint speed 


using one-way ANOVA, with sex included 


min 


to body temperature, respectively. With the best fitting, 
the modified Gaussian regression was employed to 
analyze the sprint speed, with the following equation: 
Sprint speed = y? + ae ^ {-0.5 [(T, - b)/c]}; where T, is 
body temperature; a, b, c and y? are parameters; and 
e is the natural constant. This was used to calculate 
the optimal temperature for sprint speed. Repeated 
measures ANOVAs were used to analyze the effect of test 
temperature and sex on resting metabolic rate and sprint 
speed. 

Data of previously illustrated traits of T. septentrionalis 
from Zhoushan, Zhejiang (29°32’-31°04’N, 121°30’- 
123°25’E; ~200 m), T. wolteri from Chuzhou, Anhui 
(32°15’-32°21’ N, 118°07’-118°18’ E; ~260 m), and 
T. sexlineatus from Shaoguan, Guangdong (24°50’ N, 
113°30 E; ~290 m) (Ji et al., 1996a; Chen et al., 2003; 
Zhang and Ji, 2004; unpublished data) were collected 
from published figures with Plot Digitizer (http:// 
plotdigitizer.sourceforge.net/). The thermal traits of T. 
septentrionalis (collected throughout May and July), T. 
wolteri (collected in May), and T. sexlineatus (collected in 
mid-April) were obtained from the literature. Only data 
collected during the reproductive season were selected to 
ensure comparability with the results of the current study 
(Ji et al., 1996a; Chen et al., 2003; Zhang and Ji, 2004). 
General linear models were fitted to assess the latitudinal 
tendency of Ter CT ae CT min, and T 
of Takydromus lizards. As the four species inhabit similar 


opt for sprint speed 
altitudes (i.e., from ~200 m to ~290 m), we did not include 
this variable in the model. The fact that elevation is fixed 
in the design of the study provides clarity to be able to 
detect the effect of latitude. 


3. Results 


The Tep CT nine and CT max for T. amurensis were 32.55 + 
0.46 °C, 3.06 + 0.13 °C, and 45.10 + 0.11°C, respectively 


(see details in Table 1). The sprint speed was significantly 


max 


affected by body temperature, increasing with body 
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temperature from 18 °C to around 33 °C and then 
decreasing at higher temperatures (F44 = 5.495, P < 
0.0001; Figure 2a). The optimal body temperature for 
T. amurensis that maximized sprint speed was 33.20 °C 
according to the modified Gaussian regression. Resting 
metabolic rate was significantly positively related to 
temperature from 18 °C to 38 °C (F44, = 40.708, P < 
0.0001; Figure 2b). The Sex had no effect on the resting 
metabolic rate or sprint speed (both P > 0.05). 

In the four Takydromus species, Tı were similar, as 
indicated by the non-significant tendency along latitudes 
(y = 32.99-0.55x, R’ = 0.37, P = 0.38; Figure 3a). When 
we integrated the published thermal tolerance data from 
other Takydromus lizards, we observed that the CT max 
values (y = 45.43-0.93x, R’ = 0.78, P = 0.11) remained 
similar; however, CT „in decreased as latitude increased 
(y = 2.137 + 1.04x, R’= 0.95, P = 0.02) (Figure 3b), 
indicating that the thermal tolerance range increased with 
latitude. 


Sprint speed (m/s) 


Metabolic rates (mL CO>/g/hr) 
o 
S 


18 23 28 33 38 


Test temperature (°C) 


Figure 2 Locomotor performance (sprint speed) (a) and resting 
metabolic rate (RMR) (b) of T. amurensis at different temperatures. 
RMR was expressed as CO2 production per gram body mass per 
hour (mL/g/h). Data are expressed as mean + SEM. 
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Table 1 Snout-vent length, body mass, selected body temperatures, 
the critical thermal maximum, critical thermal minimum and 
thermal tolerance range (TTR) of Takydromus amurensis. The 
sample size for all traits are equally 13. 


Mean SE Minimum Maximum 
BM (g) 3.627 0.216 2.591 4.888 
SVL (mm) 62.923 0.677 59.000 68.000 
Ty FC) 32.550 0.459 30.100 35.950 
CT ral PC) 45.102 0.109 44.510 45.845 
CT 4, lC) 3.062 0.133 2.407 4.153 
TTR (°C) 42.040 0.156 40.871 42.822 


4. Discussion 


In the context of climate change, understanding the 
thermal biology traits of animals is fundamental to 
determining organismal vulnerabilities (e.g., Sunday et 
al., 2011, 2012, 2014). In the current study, the thermal 
biology traits of T. amurensis were quantified at the 
beginning of the reproductive season, to facilitate a 
comparison with available data of congeners collected 
using similar methods. The Tep CT miwn CT maw and sprint 
speed among different species of the genus Takydromus 
were shown to be comparable and ecologically 
meaningful (Ji et al., 1995, 1996a; Du et al., 2000; Zhang 
and Ji, 2004). Since the four populations live at similar 
latitudes (i.e., ranged from ~200 m to ~290 m), we only 
tested the effects of latitude. 

T.a is the body temperature at which biological 
processes function at an optimal/suboptimal level (Van 
Damme et al., 1991; Hertz et al., 1993; Blouin-Demers 
et al., 2000; Angilletta et al., 2002a; Ortega and Martín- 
Vallejo, 2019). In T. amurensis, the optimal temperature 
(33.20 °C vs. 32.55 


°C), indicating that T. amurensis can accurately regulate 


for sprint speed was similar to Ta 
its body temperature to perform optimal/suboptimal 
functions. However, as we did not measure operative 
temperature (T,) and active body temperatures 
(T,), we were unable to determine the efficiency of 
thermoregulation. Numerous reptiles have been reported 
to have temporal, spatial, and individual variation in T,,, 
(Gatten Jr, 1974; Andrews, 1998; Stellatelli et al., 2018; 
Refsnider et al., 2019). There is also increasing evidence 
demonstrating that thermal biology traits vary among 
seasons (e.g., Yang et al., 2008; Sun et al., 2014; see details 
in below). However, Tsel is normally only determined 
during the reproductive season and assumed representing 
the Tsel of the given species or population (Ji et al., 1995, 
1996a; Shine and Madsen, 1996; Corovié and Crnobrnja- 
Isailović, 2018; Refsnider et al., 2019). It should be noted 


Selected body temperature (°C) 


Body temperature (°C) 


T. wolteri 


T. amurensis 


T. septentrionalis T. sexlineatus 


Figure 3 Tı (a) and Cimas and Cymin (b) of four Takydromus liz- 
ards. Data are expressed as mean + SEM. Species are listed from 
high to low latitude in X axis. Data are collected by Plot Digitizer 
(http://plotdigitizer.sourceforge.net/), from the figures and avail- 
able data of references (Ji et al., 1996a; Xu and Ji, 2006; Zhang 
and Ji, 2004). 


sel Tmin 


that although T. amurensis is a cold-climate species, the 
T,«a in T. amurensis is very similar to published data for 
T. wolteri, T. septentrionalis, and T. sexlineatus (Figure 3a). 
Additionally, the lizards Zootoca vivipara, Lacerta agilis, and 
Iberolacerta bonnali which also inhabit cold climates, use 
thermoregulation to maintain their body temperature 
warm in responding to cold climates (Herczeg et al., 2003; 
Yang et al., 2015; Ortega et al., 2016). Influence of the 
thermal environment and activity patterns could induce 
variations in T,,, (Ji et al., 1996a, 1996b; Du et al., 2000; 
Zhang and Ji, 2004). For example, lizards in forested 
habitats are likely to have lower T,,, than those using 
more open habitats. The T,., of Sphenomorphus indicus, 
whose habitats are forested, was as low as 25.7°C (Ji et 
al., 1996b). On the contrary, the lizards using more open 
habitats and employing basking behaviors would exhibit 
higher selected body temperatures (e.g., Li et al., 2017; 
Wang et al., 2019). T. amurensis is usually active in open 
microhabitat, with shuttling between open and shaded 
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patches for thermoregulation (Zhao, 2002), which may 
plausibly explain high T,« of T. amurensis. Admittedly, we 
conducted selected body temperatures measurements 
of T. amurensis with established method according to 
those used in other Takydromus species (e.g., Ji et al., 1995; 
Zhang and Ji, 2004), increasing literatures indicate that 
the selected body temperatures of ectotherms should be 
determined with multiple repeated measurenets (i.e., eight 
to ten measurements per individual) (e.g., Hertz et al., 
1993; Li et al., 2017). We also encourage further researches 
should employ multiple measurements to make a better 
understanding of thermal preferences in lizards, including 
Takydromu lizards. 

The thermal tolerance range increased toward high 
latitudes in Takydromus lizards, with relative conservative 
CT max and decreased CT min towards high latitudes. (Figure 
3b). Microhabitat usage can significantly influence CT 
and CT min: For example, CT max was lower in S. indicus, 
which uses shaded habitats likely to be cooler (37.6 °C; 
Ji et al., 1996b) than species occupying open habitats, 
such as Eumeces elegans (41.9 °C; Du et al., 2000) and T. 
septentrionalis (42.3 °C; Ji et al., 1996a). Among sympatric 
lizards, Eremias argus (45.8 °C) and E. multiocellta (45.1 °C) 
that occupying filtered microhabitats, had lower CT max 


max 


min’ 


than Phrynocephalus przewalskii (47.1 °C), which uses open 
habitats (Li et al., 2017; Wang et al., 2019). Analogously, 
CT nin Was noticeably greater in lizards that use open 
habitats. The CT „in of E. chinensis (6.3 °C; Ji et al., 1995) and 
E. elegans (9.3 °C; Du et al., 2000), species occupying open 


min 


habitats, were greater than species using more spatially 
complex habitats, such as E. argus (1.0 °C; Luo et al., 2006). 
Although all the four species compared in the current 
study use open habitats (Zhao et al., 1999), T. amurensis 
may spend more time on open habitats in response to 
lower temperature conditions, as happens in other high- 
latitude distributed lizards (Sunday et al., 2011; Sunday 
et al., 2012). This activity pattern may result in increased 
time invested in thermoregulation (i.e., basking) and thus 
result in higher CT max Based on this assumption, global 
meta-analyses have demonstrated that CT „a increases 
with latitude in lizards (Sunday et al., 2011; Sunday et 
al., 2012). In contrast, CT „in is lower at high latitudes, 
likely because species inhabiting lower latitudinal or 
warmer climatic regions do not face the challenges of 
cold and seasonally varying temperatures (Janzen, 1967; 
Ghalambor et al., 2006). Although we studied the thermal 
biology traits of Takydromus lizards across a wide latitude 
range, our results are preliminary and can be improved 
upon by incorporating more species and populations. 

The resting metabolic rates of T. amurensis were 
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thermally dependent, with allometric enhancement as the 
temperature increased from 18 °C to 38 °C (Figure 2b). 
As the basic biological rate, the resting metabolic rates of 
lizards vary among species, populations, and individuals. 
The resting metabolic rates of T. amurensis were lower 
than those of the desert-dwelling lizards E. argus, E. 
multiocellata, and P. przewalskii at even temperatures, 
although they are distributed at similar latitudes (Li et 
al., 2017; Wang et al., 2019). Unfortunately, the resting 
metabolic rate of T. septentrionalis, T. wolteri, and T. 
sexlineatus has not been studied extensively; this means 
we were unable to make comparisons among Takydromus 
lizards. The sprint speeds of T. amurensis and other 
Takydromus species were enhanced as the temperatures 
increased until the optimal temperature was reached, 
and then decreased as the temperature increased beyond 
the optimal temperature (Ji et al., 1996a; Zhang and 
Ji, 2004). However, in all four species we compared, T. 
amurensis had the lowest sprint speed (Figure 4a). Notably, 
the optimal body temperature for the sprint speed of T. 
amurensis was 33.20 °C. In the other three species, the 
optimal temperatures for sprint speed increased as the 
latitude decreased (Figure 4b). For species from warm 
and stable thermal environments, sprint speed is more 
thermally dependent, and the optimal temperatures are 
higher, as shown in T. sexlineatus (Figure 4; Zhang and Ji, 
2004). Why T. amurensis, a cold-climate species, has such a 
high optimal temperature for sprint speed is still unclear. 
However, the optimal temperature for sprint speed (i.e., 
33.20 °C), which is roughly equivalent to T,,, (i.e. 32.55 
°C), imposed an optimal/suboptimal sprint speed at T,« in 
T. amurensis. Future studies could investigate the effect of 
temperature at a finer scale (e.g., 2 °C between each test 
temperature), to further clarify the thermal dependence 
of sprint speed in T. amurensis and thus determine the 
optimal temperature more precisely. Although the 


mechanism is unclear, Tı and To might be subject to 


opt 
both thermal environment and confounding evolutionary 
factors (Huey and Kingsolver, 1989). Alternatively, 
maintaining T, at a level that appears suboptimal may 
allow the animal to remain responsive to environmental 
variation and, if required, increase activity such as those 
associated with foraging and escaping, following the 
‘suboptimal is optimal’ hypothesis (Martin and Huey, 
2008). For example, the range of T, and T, for some 
Liolaemidae lizards is below both their preferred body 
temperature and the optimal temperature at which they 
can reach maximum locomotor performance (Bonino 
et al., 2011). Addtionally, T, at suboptimal temperatures 
is also beneficial to avoid over-heating effects in 
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performance, since the performance would drop sharply 
if the optimal temperatures are surpassed (e.g., Martin 
and Huey, 2008). 

Based on our current preliminary pattern of T 


CT 


is plausible that T. amurensis would be less vulnerable to 


CT 


and T» for sprint speed in Takydromus lizards, it 


sel? min’ 


max? 


the warming temperatures imposed by climate change 
than other Takydromus species in habitat in low latitudes. 
First, the CT nax 


indicate that the thermal safety margin is large. However, 


for T. amurensis is high; thus, this may 


we were unable to verify this as we did not calculate the 
thermal safety margin owing to the lack of precise data on 
environmental temperatures (Sunday et al., 2014). Second, 
as T. amurensis inhabit areas with cooler temperatures 
than their physiological optima, it is possible that 
warming temperatures may enhance their fitness, as in 
high-latitude species (e.g., Cabezas Cartes et al., 2019; 
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Figure 4 Sprint speed of Takydromus lizards (a) and their Topt (b). 
(a) Arrows indicate the Topt for sprint speed, respectively. Data 
are expressed as mean + SEM. Species are listed from high to low 
latitude in X axis. Data are collected by Plot Digitizer (http:// 
plotdigitizer.sourceforge.net/), from the figures of references (Ji et 
al., 1996a; Xu and Ji., 2006; Zhang and Ji, 2004). 
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Gomez Alés et al., 2019). 

The thermal biology traits of T. amurensis were collected 
during the reproductive season (i.e., late April). However, 
it is possible that seasonal plasticity or acclimatization/ 
acclimation can influence thermal biology traits and their 
associated performance curve in reptiles and thus cause 
variations of thermal biology traits in Takydromus lizards. 
For example, Takydromus septentrionalis from different 
locations did not differ in their selected body temperatures 
and thermal tolerance after identical thermal acclimation 
(Yang et al., 2008). Furthermore, the thermal performance 
curve of the sprint speed of Plestiodon chinensis was shown 
to vary between seasons (Sun et al., 2014). The metabolic 
enzyme activity in the muscles of Alligator mississippiensis 
also responds to temperatures differently between 
seasonal acclimation (Seebacher et al., 2003). Determining 
the seasonal variation of biological traits in cold-climate 
species such as T. amurensis is important because thermal 
fluctuations are predicted to be more significant at high 
latitudes (IPCC 2013). Another potential cause of the 
observed variation in the thermal biology traits of T. 
amurensis in the current study may be the limited sample 
size (i.e., n = 13). Larger sample sizes of individuals will be 
required to investigate seasonal acclimation and confirm 
our preliminary results on latitudinal trends. However, 
as a cold-climate species, the natural population density 
of T. amurensis is low and is threatened by habitat 
destruction (Zhao, 2002; Portniagina et al., 2019). In 
future studies, meta-analysis methods may be more 
suitable for revealing the geographical patterns of thermal 
biology traits within congeners. It is possible that the 
inter-specific variation in thermal biology traits along 
latitudinal gradients is attributed to other factors, such as 
precipitation, vegetation cover, or food availability, which 
also varies across altitudes. However, as the Takydromus 
species occupy similar elevations, it is plausible that the 
observed pattern among species is induced by latitudinal 
differences. Future studies, investigating more species 
distributed over a large geographical span, should be 
conducted to test altitudinal pattern. 

Our study determined the thermal biology traits of 
cold-climate distributed T. amurensis, summarized the 
interspecies variation, and provides preliminary findings 
suggesting that T CT mav and CT, 
gradients in Takydromus lizards. These findings will 


min Vary along latitudinal 
contribute to the understanding of the thermal adaptation 
of reptiles and establish basic criteria for determining 
the climatic vulnerability of Takydromus species. Future 
studies should investigate more details of the biological 
response to thermal environments, including acclimation 
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and life history variations. Furthermore, building on the 
data collected in this study with additional measurements 
of operative temperature and active body temperature 
in field, we will be able to calculate the thermal safety 
margin and the effectiveness of thermoregulation, which 
would be helpful in evaluating the degree of safety or 
beharioral flexibility to thermal variation, in the context 
of climate warming. Additional research into thermal 
biological responses will facilitate projection of the 
vulnerabilities of Takydromus species to climate warming. 
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